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1. INTRODUCTION 

The realization of two-component condensates jl| as 
well as condensates in a double-well potential [pjhas 
attracted considerable interest, both experimental and 
theoretical (see, e.g., [|J and references therein) to the 
phenomena caused by phase coherence of two conden- 
sate modes. In the work Q spatial quantum coherence 
was observed by means of interference pattern in two 
overlapping condensates. This interference pattern was 
conformed in the work Q with the numerical simulation 
of the Gross-Pitaevskii equation. In the work || coher- 
ent oscillations of the relative populations have been ob- 
served in driven two-component condensates with differ- 
ent internal states. As is well known, a bright manifes- 
tation of phase coherence is Josephson effect. In numer- 
ous studies devoted to the Josephson effect in systems 
of two condensates in different internal states ||, || or 
in a double- well potential j7[-|| coherent Josephson os- 
cillations are considered for various dynamical regimes 
caused by the competition between tunneling and intra- 
condensate interaction (nonlinearity) . In the work 
the Josephson coupling energy is calculated for small- 
amplitude oscillations in a double-well potential. There 
are estimated damping effects due to the normal cur- 
rents at finite temperature. In the work || it is shown 
that for a relatively weak interaction, the particle num- 
ber oscillations between the condensates are complete. 
They are suppressed when the total number of atoms in 
the condensates exceeds a critical value and the behavior 
of the system is governed by the nonlinearity. Nonlin- 
ear Josephson-type oscillations in the relative oscillations 
of driven two-component condensates are studied in the 
work [g] . Decoherence effects and quantum corrections to 
mean- field solutions have been considered in ||] , |l(| . In 
the damping effects (even at zero temperature) of the 
Josephson current are derived by means of an functional 
integral approach. The detailed treatment of the nonlin- 
ear classical dynamic of the condensates in a double- well 
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potential has been performed in the works p2[,p^|. In 
the work fl4|| the quantum and thermal fluctuations of 
the phase are studied for the condensates in the double 
well potential. 

It should be emphasized that experimental observation 
of the Josephson effect is difficult because the small en- 
ergy splitting associated with Josephson coupling means 
that thermal and quantum fluctuations will destroy the 
phase coherence between two condensates even at the 
lowest achievable temperatures |ll|, jl4|]. While the en- 
ergy splitting can be increased, for instance, by lowering 
of the barrier height, it then becomes comparable with 
that of motion states of the condensates. 

However, the problem of the interaction between the 
Josephson degrees of freedom and states of motion (oscil- 
lations) of the trapped condensate has yet to be analyzed. 
The present paper focuses on a mechanism of formation 
of the bound state of the Josephson degrees of freedom 
and trap oscillations due to their interaction. The mech- 
anism proposed in what follows may be important for 
detection of the excited Josephson states. 

The results achieved in this paper may be listed as 
follows. 

1. To formulate the problem of the interaction between 
the Josephson and oscillation degrees of freedom in an ad- 
equate manner, the quantizied spectrum of the particle 
number generated by Josephson coupling is derived. The 
states of this spectrum represent the quantum analogue 
of the nonlinear coherent Josephson oscillations consid- 
ered in p2| . In what follows, the states of the quantizied 
spectrum are named as Josephson states. The spectrum 
is highly nonequidistant and has the logarithmic singu- 
larity in the density of states at the 2Ej energy, where 
Ej is the Josephson coupling energy. As is shown in what 
follows, any Josephson state can be realized by means of 
a given initial disbalance of the particle number in two 
condensates. 

2. There is considered the interaction between the trap 
oscillations and excited Josephson states corresponding 
to the large enough initial disbalance of the particle num- 
ber. It is shown that this interaction is responsible for 
the formation of the bound state of n m ^> 1 oscillation 
quanta with the Josephson state corresponding to the 
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definite initial disbalance. In the Thomas-Fermi approx- 
imation at > w , where fj, is the chemical potential 
of the condensates, ujq is the characteristic frequency of 
a trap, h = 1, the bound state arises in the region of 
the dense oscillation spectrum. In this region the level 
separations are small in comparison with the harmonic 
oscillation frequency loq. The equilibrium values of the 
oscillation quanta and initial disbalance of the particle 
number are coupled self-consistently and can realize an 
absolute minimum of the thermodynamic energy at the 
large enough interaction. The thermodynamic average 
n m 7^ generates the equilibrium distortion of the con- 
densate shape. It allows us to detect the Josephson states 
by observing the change in the condensate shape. 



THE QUANTUM SPECTRUM OF THE 
PARTICLE NUMBER 



A Josephson coupling is realized for the condensates 
in an symmetrical/asymmetrical double- well potential 
formed by two different traps with a barrier between 
them The barrier is created by the laser light, and 
its height is directly proportional to the laser power and, 
thus, can be varied with ease. The proposed mech- 
anism is also suitable for the condensates in different 
internal states in the same trap. Experimentally, this 
may be a superposition of two Rb 87 condensates in the 
\F = l,m F = — 1>, j-F 7 = 2,m F = 1} states §, [§. In 
the case of a weak Josephson coupling, the basis states 
are the self-consistent ground states in two condensates 
separately. The wave function of the condensate with 
Josephson coupling takes the form of a superposition of 
these states, namely, ^(r,t) = if>i(r)ai(t) + ^2(r)a 2 (t), 
where ipi(r) are the normalized solutions of the Gross- 
Pitaevskii equation, ai {t) = rf /2 (t)e ie *W, i = 1,2, 2V 1)2 , 
9i (t) are the particle numbers and phases of each conden- 
sate. 

As is well known ||, @], §, (l3) the Hamiltonian of 
two condensates with Josephson coupling has the form 



Hj -E = E C (AN) 2 - 2Ej cos < 



(1) 



here AJV = Ni - N 2 , E c = dfi/dN, fj, = m = p, 2 are the 
chemical potentials, N are the total particle numbers, 
4> = Oi — 9-2 is the relative phase of the condensates. The 
quantities Ec, Ej depend on the total particle number 
N . (AN), 4> are canonically conjugate variables. In (|l|) 
the energy origin is the mean-field summary energy of 
the condensates, namely, Eq = fiN. 

The quantization of the Hamiltonian (0) produces the 
spectrum of the particle number in the Josephson poten- 
tial Ej cos <j). As is shown in what follows, any Josephson 
state can be realized by means of a given initial disbal- 
ance of the particle number. For this reason, it is inter- 
esting to obtain the complete spectrum generated by the 
Hamiltonian of ([!]), and to represent it as a function of 
the initial disbalance. 



The Schrodinger equation for the Hamiltonian ([l]) is 
derived by means of the quantization rule: (AN) — > 
—id/d<f). As a result, we obtain Mathieu's equation 



2Ej cos < 



* = e*, e = Hj - E , (2) 



For e > 2Ej this equation has a continuous spectrum. 
The states of this spectrum correspond to the classical 
states with unlimited phase change: — oo < 4> < +oo that 
are named as the self-trapping states in the works JlJ, [l3) . 
In the region —2Ej < e < 2Ej Eq.(||) has the discrete 
spectrum. It corresponds to the region of finite motion of 
Hamiltonian (Q), where the relative phase changes within 
— arccos(e/2£'j) < <p < &rccos(e/2Ej) for every e. In the 
Josephson regime at Ej » Ec || , the number of levels 
in a well is large and the discrete spectrum is specified 
by the Bohr-Sommerfeld formula: 



v(e v ) = & -^-AN(cp r ;e„) 



— (e v + 2Ej cos < 
Ec 



v c E(k) - (1 - k 2 )K(k) ■ 



1/2 



(3) 



8(Ej\ 1/2 2 e + 2Ej 

K(k) and E(n) are complete elliptic integrals of the first 
and second kinds. The density of states pd(s v ) results 
from Eq.(|5|) and is equal to 



1 dv K(k) 
Pd(e„) = 7T1~ = ~ 



(4) 



The level separations in (Q) are given by the following 
expression 

"<«•>-§- - = 2 <W 2 - < 5 > 

At (Ej/E c ) 1/2 > 1 we have the relation uj m <C 2Ej. 
The quantity oj m determines the maximum splitting of 
the levels in the Josephson well. In what follows, the 
states with v < v c are denoted as "libration" states. 

In the region of v > v c the v(e) dependence and density 
of states are determined by the following expressions 

TT Z \E C J 2tt z UJ m 

_ (6) 
The states (g) with v > v c are named as the self-trapping 
states. 

Eqs.(|), (|) imply that 



e(v) w ~2Ej + uj m v; 1 < v < v c , 
e(v) w Att 2 E c v 2 ; V»v c . 



(7) 
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At the same time, it is easy to show that d 2 e(v) / dv 2 < 
at v < v c and d 2 e(y)jdv 2 > at v > v c . At v = v c the 
curve e(v) has a flex point. 

Since the energy is conserved, the state with a given 
v value can be realized by means of the definition of 
the initial values of (A7V)o and </>(0). Namely, e{y) — 
E c (AN)l - 2Ejcos(f>(Q), Supposing that 0(0) = 0, we 
obtain the following relation between v and (AiV) 

e{y) = -2Ej + E c (AN)l (8) 

Using Eqs.(0), (j|), we arrive at the expressions 



v = — — | (AiV)o | ; for the self-trapping (sf) states, 
— ) (AN)l; for the "libration" (1) states. 



(9) 



Combining (Q), (||), we come to the depedence 

p d , c (e) ~u m 1 \n\l-e/2Ej\- 1 ; e -» 2E± (10) 



Thus, there appears a new logarithmic singularity at the 
boundary separating the libration and self-trapping spec- 
tra. 



3. THE INTERACTION OF THE JOSEPHSON 
STATES AND CONDENSATE OSCILLATIONS 

In this and next sections it is shown that the spectrum 
of the system can change drastically due to the interac- 
tion between the excited Josephson states (Q)-(||) with 
large enough v values and oscillations of the condensate. 

The interaction can be realized by the following mech- 
anisms. First, the interaction is secured if we allow for 
the dependence of Ec in Eq.(||) on the atom displace- 
ments. The latter are generated by the condensate oscil- 
lation. Second, the interaction can be realized by apply- 
ing two-photon traveling-wave laser pulse with the Rabi 
frequency fl. The pulse both creates the condensates 
with the different particle numbers and induces the inter- 
action of atom's displacements with the excited Joseph- 
son states corresponding to the particle number disbal- 
ance created by the pulse. General description proposed 
in what follows is independent on the specific mechanism 
producing the interaction. 

Let us consider the classical states of motion of the 
condensate. These states may be described in terms 
of the complex amplitudes a*, a = n 1 / 2 e ±I¥ ' 1 , where 
n = (a\a + a\a) = \a\ 2 is the average number of quanta in 
the coherent state a). The variables n, (p% are canonical. 
By classical state of motion we mean that it's number of 
quanta is very large, n 3> 1. It is convenient to specify 
the relation between the a, a* amplitudes and a, a + op- 
erators by the following way: a = N^ 1 / 2 ^. At this, the 
commutator of a, a* is equal to zero with macroscopic 
accuracy: [a, a*} = 1/N — > 0. The Hamiltonian of the 
motion states can be written in the form: Ne(n). 



For the quasiclassical Josephson states with v ^> 1, 
the c„ , c* amplitudes may be written in the form c„ — 
v l / 2 e %tp<2 . However, it is convenient to rewrite e(y), c„ in 
terms of the variable x = |(AiV)o|/./V 1/2 > 1. Com- 
bining this inequality with the requirement that x — 
\{AN)o\/N « 1 we arrive at the conditions for the x 
values: 



1 < x <C N 1/2 . 



(11) 



Using Eqs.(0),(||), we find that v — v(x), e(v) — 
N(— Ej /N+Ecx 2 ). In general case.the v(x) dependence 
is implicit. It is determined by Eqs.(||), (g), (|j). However, 
in the particular cases of the "libration" (e(v) <C Ej) and 
self-trapping (e(v) 3> Ej) states, the relations between 
v and (AAQp can be represented in a simple form, as is 
seen from (g). By means of Eqs.(||), @, we come to the 
following expressions 



H Q = Ne {n, x) = N[--± + e(n) + E c x 2 \ 



N 



= N[-^+e(n)+ej(x)], 



(12) 



M =N (E£.) 1/2 xe i<p> = N (E£ 



1/2 



(13) 



For any mechanism producing the interaction between 
two subsystems, it can be written in the form of a mul- 
tiple Fourier series in tpi, (f2 



H int = NY, l9k lk2 (N)a* k 'c k u 2 + c.c] 

ki k,2 



+ c.c. 



(14) 



fci ,k 2 



Here a s f — 1/2, ai — 1, k\,k2 are integer. For the sake 
of simplicity, we disregard the phase-independent inter- 
action. Using Eq.(n3h, we obtain 



fc 2 /2 



The constant g is specified by the concrete mechanism 
producing the interaction. Let us assume that the term 
with the phase <jf k — k\ r tp\ — k2r<P2, which varies anoma- 
lously slowly with time, can be set off in sum (|l4|). It is 
possible to make under two conditions 



, ,' M n ) \ ( dej(x) \ 

Kir — i = k 2r I ; I , ^ 



dn 



dx 



de(n) 
dn 



dx 2 



(Ax) 

max 



2k 2r E c ' 

( dHi n t 



\ 9X Jx 



(16) 



(17) 
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Eq.(p7|) is written with taking into account that H , H int 
are functions of a single dynamic variable, e.g., x. The 
e' n quantity defines the level separations of an oscillation 
spectrum. 



As is shown in what follows, the condition (16) is equiv- 
alent to that of the minimum of the function £o(n, x) over 
x. In turn, when the minimum exists, it can provide the 
principal contribution into the thermodynamic functions. 

In addition, the condition ( |l6| ) implies that the phase 
fa is an approximate integral of motion in the ab- 
sence of changing x near x m : d<ff k /dt w k lr dH / dn — 
k 2r dH /dx w 0. 

The inequality (|l7|) implies that the width of the near- 
minimum region is large at the characteristic interac- 
tion variation scale. From Eqs.(fl6l), (|l7| ) one can ob- 
tain that time changing the fa phase is proportional to 
(d 2 Ho/dx 2 ) m Ax, where Ax is the change in x near the 
x m value. The maximum value (Ax) max specifies the 
width of the near-minimum region in that dfa/dt ~ Ax. 
The estimation for (Ax) max is given in what follows. 
Thus, the leading term in sum ( |l4| ) has the form 

H int = N 9k' l \ N ' n > X ) C °S^fc, fa = felr^l - k 2 r^2] 

All remaining terms in this sum are rapidly oscillating 
perturbations and will be disregarded in this work. Here 
and below the index k in g k and fa denotes a set of 
k\ r ,k,2r- One can easy to show that aside from the 
energy H = N[eo(n,x) + g k (N; n, x) cos fa)] the sys- 
tem in question has an additional integral of motion: 
no = n/kvr + x/k2n dng/dt — 0. Owing to this, the 



condition (16) is equivalent to that of the minimum of 
eo(no,x) over x at a given no value, as is mentioned 
above. Using Eqs.(|ilf)-(|l8|), it is straightforward to write 

(r) 

the Hamiltonian H m = Hq + A near the minimum to 
the first nonvanishing order in Ax: 



H m = N 



e (n m ; x m ) + 



d 2 e 
dx 2 



(Ax) 2 



fjkr 



. cos fa 
(19) 



n rn = n - k lr x m /k 2r , (d 2 e Q /dx 2 ) m = 2E C , g km = 
gjff' l '(N;n m ,x m ). Terms with the derivatives of Hi nt 
are absent in Eq.(|l9|) due to the condition (|l7|). 

Using that E c ~ uWa/a ) 2 / 5 A^ 3 / 5 in the Thomas- 
Fermi approximation |L5| (here a, ao are the scattering 
and oscillator lengths, respectively) , we can represent the 
range of 1 <C x m <C N 1 / 2 in the form 



1 (Na 




(20) 



As is known Jig] , the relation (Na/ao) 1 oc- 

curs in the Thomas-Fermi approximation. However, 
N^ 1 / 2 (Na/a ) 2/r ° <C 1. Owing to this relation, the con- 
dition ( |20| ) (or, what is the same, the condition (|l6|)) 
specifies the region of the dense oscillation spectrum, 
where e n <C wo- Both here and in what follows, we sup- 
pose that k\ r — ki r — 1 for the sake of simplicity. 



From Eq.([T^) the (Ax) max can be estimate to be 
(Ax) max ~ {gkm/E c ) 1/2 - Hence, the condition (0) 
takes the form: 



E c 



gkm(N) 



1/2 



(21) 



In what follows (see Eq.([34|)), it will be shown that the 
relation (Elh is fulfilled with macroscopic accuracy. 



4. THE GROUND STATE 

At a fixed uq value the principal contribution to the 
partition function comes from the neighbourhood of the 
minimum at x = x rn . The expression for Z(hq; x rn ; T) is 
equal to 



OC 7T 

Z (no; x m ;T) = const J d(Ax) J dfat 



-/3H m (n ,(Ax),<j> k ) 



const 



= / /Q at r~> m/2 ex P [-/3Ne (N; n m ; x m ) + In I (/3Ng km )] , 
(fJNE c ) y l l 

(22) 

where f3 = 1/T, T is temperature, Iq(x) is the modified 
Bessel function. Eq.(|2^) implies that the expression for 
the free energy of the system has the form 

F = Ne (n m ;x m ) + ^Tln(f3NE c ) - T In I Q ((3N g > km ), 

_ " (23) 

Using Eq.(|2^), we come to the following equation for 
the fi m value realizing the minimum of the free energy 



dear. 
dn„ 



dg k , 
dn„ 



h(/3Ng km ) 
Io((3N9km) 



, (24) 



here I\(x) — Iq(x). In addition to n m , the thermodv 
namic average of cos fa may be determined from Eqs.(|22 
or (p3|) . This average is equal to 



(COS fa)rj 



d\nZ 



dF 



h(f3Ng km ) 



d(fiNg km ) Td((3Ng km ) I ((3Ng km )' 

(25) 

The order parameters n m , (cos fa) t describe new coher- 
ent state. There is the bound state of the n m oscillation 
quanta and Josephson state generated by the initial dis- 
balance of the particle number that corresponds to the 
x m value. In addition, this state has the equilibrium 
phase coherence factor (cos fa)T- The fi m ^ value pro- 
vides the equilibrium distortion of the condensate shape. 
The aboveobtained equations imply that the shape dis- 
tortion is sclf-consistcntly coupled to the x m value defin- 
ing the equilibrium initial disbalance of the particle num- 
ber. 

At T = the n m value realizes the minimum of the 
thermodynamic energy 



E,, 



N[e (N; )]■ (26) 
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To determine n m (T = 0), it is suitable to use the follow- 
ing consideration. As is well known, the level separations 
e' n change slowly in dependence on n within the range of 
dense (quasiclassical) spectrum. For this reason, we can 
suppose that e n « const = LUb <C ojq. On this assump- 
tion, the x m value does not depend on n m and the fi m 
quantity is equal to 



,1/2 _ 9 



(sf,l) 



1 



a sf,l 



(27) 



2uj h ' 2' ' 

(It is worth noting that (d 2 E m / 'dn^J > 0). Both 
here and in what follows, the denotation (N) = 

g( s f' l )(N) is used. Taking Eq.(^7|) into account, one read- 
ily gets 



=(•/) 
9lm 



fJ 



(s/)2 



2wb 



fflm 



,(02 



(AO* 2 , 



2^6 



The transition to the state with n m ^ 0, (cos</)fc)T 
is second-order. By imposing that n m — > at the 
transition temperature, one gets from Eq. (|2~i|) . 



■N. 



(32) 



The dependences of the transition temperatures ( J32] ) on 
the total particle number are given by the expressions 



(33) 



Thus, the transition temperature T c has the macro- 
scopic smallness in comparison with the T c tempera- 
ture. Along with the conditions (|30|), (|3l|), this fact im- 
plies that the libration Josephson state forms the bound 
state with the condensate oscillation rather than the self- 
trapping state. 

5. CONCLUDING REMARKS 



E, 



N 



E, 



(0 



N 



_Ej_ 
N 

El 

N 



9 



(sf)2 



(AT) 



fJ 



8Ec 
^ 2 {N)x r , 



16E C 



*E 2 x 2 m \ 

gW) 2 (N) J 



l6E 2 x rn 



ff(')2(JV) 



(28) 



The expressions for energies imply, first, that we obtain 
the minimum in the region of the dense en oug h oscillation 
spectrum, which satisfies the condition (|2y) . The min- 
imum corresponds to the formation of the bound state 
for the n mi x m values. Second, as is seen from Eq.(28), 
the absolute minimum of E m can be realized within the 
ranges of 



,(*/) 



1 < x m < 



(AO 



8E C 



1 < x m < 



9 {l)2 (N) 
16E 2 ■ 



(29) 



These conditions are met when the interaction matrix 
elaments g^f'^^N) are large enough. Let us estimate the 
condensate parameters that are required for existence of 
the absolute minimum. In the Thomas-Fermi approxi- 
mation, the inequality g^ s ^(N) » Eq occurs provided 
the total particle number is not very large, namely, 



9 /ao\ 8 / 3 
N < — — 
u>o V a ) 



(30) 



In turn, the relation g^ 2 (N) 3> Eq holds true within 
the range of 



fi 1 /2 w 3/2 



V a 



0.6 



»1. 



(31) 



Here we use that Ej = Q.N. The condition (|3l]) is ful- 



filled for all admissible parameters, if g 2 /(£l 



1/2,3/2 



It should be emphasized that the right-side inequali- 
ties in Eq.(p9|) are much stronger than the condition 

x m «A^ 1 /2. 



We have found that the interaction between the 
Josephson and oscillation states results in new coher- 
ent ground state of the double-well trapped condensate. 
There is the bound state of the anharmonic condensate 
oscillation and excited Josephson state. The latter is 
specified by the definite initial disbalance of the particle 
number. The x m value defining the disbalance is self- 
consistently coupled to the equilibrium number of the 
oscillation quanta entering the bound state. Along with 
the n m order parameter, this state has the equilibrium 
phase coherence factor (cos 4>k)T, where 4>k is the relative 
phase of the Josephson and oscillation modes. 

It should be emphasized that the bound state arises 
in the neighbourhood of the minimum of eo(no;x) over 
x, where the conditions (|l6|), ( p"7| ) are met. The con- 
dition ( |l6| ) specifies the type of the interaction between 
the Josephson and oscillation degrees of freedom. In ad- 
dition, it imposes the definite restrictions on the spec- 
trum of the oscillation states which may effectively in- 
teract with the Josephson degrees of freedom. Namely, 
the density of oscillation states should be large enough 
in order to satisfy the conditions (|2p|). For instance, in 
an asymmetrical double-well potential the energy e(n) of 
it's classical oscillation states has three branches. Two 
branches have energies of ei,2{n) < Vb, here Vb is the bar- 
rier height. There are the maximums at t\,i(n ma x) = Vb- 
The third branch has both the energy 63(71) > Vb and 
the minimum at e3(n m i n ) = Vb- Thus, the regions of 
the dense spectrum exist in the neighbourhood of the ex- 
tremums of the functions ei(n) (i = 1, 2, 3). In the other 
words, in the neighbourhood of barrier top. 

As is shown in section III, the condition ( |l7| ) can be 
represent in the form (21). Substituting the aboveob- 
tained expressions for fi m , g\ m into (pl[), we find that it 
takes the form 

1/2 

|(s/,0 ^ I 9\r, 



E c 



(34) 
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This inequality is satisfied with macroscopic accuracy. 

It is worth noting that the conditions (|3C|), ( |3l| ) of the 
absolute minimum of the thermodynamic energy are ex- 
perimentally controlled by means of either the interaction 
matrix element g or the particle number N. 

The existence of the bound state generates the equilib- 
rium distortion of the condensate shape specified by the 
n m value. This mechanism can provide the experimen- 
tal detection of the excited Josephson states. The latter 
can be observed by changing the condensate shape. In 
addition, the phase transition to the new ground state 



occurs at T — T c , where transition temperatures T c 
are defined by Eq. (|32|) . 
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